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Since gluons in QCD are interacting fundamental constituents just as quarks are, 
we expect that in addition to mesons made from a quark and an antiquark, there 
should also be glueballs and hybrids (bound states of quarks, antiquarks and gluons). 
In general, these states would mix strongly with the conventional qq mesons. However, 
they can also have exotic quantum numbers inaccessible to qq mesons. Confirmation 



of such states would give information on the role of "dynamical" color in low energy 
QCD. In the quenched approximation we present a lattice calculation of the masses 
of mesons with exotic quantum numbers. These hybrid mesons can mix with four 
quark {qqqq) states. The quenched approximation partially suppresses this mixing. 
Nonetheless, our hybrid interpolating fields also couple to four quark states. Using 
a four quark source operator, wc demonstrate this mixing for the 1 meson. Using 
the conventional Wilson quark action, we calculate both at reasonably light quark 
masses, intending to extrapolate to small quark mass, and near the charmed quark 
mass, where we calculate the masses of some ccg hybrid mesons. The hybrid meson 
masses are large — over 4 GeV for charmonium and more than twice the vector meson 
mass at our smallest quark mass, which is near the strange quark mass. 
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1 Introduction 



While there is a long history of glueball mass calculations in lattice QCD, including attempts 
to use lattice calculations to identify experimentally observed mesons with glueballs[0, 0, 
hybrid mesons have received much less attention. These bound states of quarks and gluons 
have been treated in a variety of approximations to QCD, including the bag model, flux-tube 
model, and QCD sum rules0. As with glueball candidates, a hybrid characterization of an 
observed state is more convincing if there is not only a match in mass, but also in decay 
branching ratios. Since the occurrence of hybrid states is obviously a nonperturbative phe- 
nomenon, and since lattice gauge theory provides an ab initio nonperturbative formulation 
of QCD, in principle lattice gauge theory is the ideal method for calculating their properties. 
One hopes that lattice methods will eventually provide reliable masses and branching ratios, 
as well as providing a basis for testing the various approximations to QCD. 

Although in principle the lattice approach is ideal, in practice there are some difficulties. 
On the lattices used here, the hybrid masses are large compared to the lattice spacing, 
so lattice spacing errors are a serious problem. Also, the propagators are noisy, although 
perhaps not so much as glueball propagators. This means that we do not have long plateaus 
in the effective mass, and we must extract mass estimates from rather short distances. While 
we will end by making the best mass estimates we can, at this stage we are still exploring 
methods and the dependence of the results on parameters such as the lattice spacing and 
quark mass. 

The earliest lattice calculations of hybrid mesons used static quarks, where hybrid states 
appear as excitations of the gluonic string[^. Also in an early study, the UKQCD collabora- 
tion studied hybrid states in the Upsilon system, in a simulation of nonrelativistic QCD[0. 
More recently, the UKQCD group has presented results in the quenched approximation for 
quark masses about equal to the strange quark mass|^, |^ and we have presented preliminary 
results using Wilson valence quarks and Kogut-Susskind sea quarks]^. 

Hybrid mesons can have the same quantum numbers as conventional qq mesons, and 
would be expected to mix strongly with them. (This mixing was demonstrated in Ref. 
In addition, hybrid mesons can have exotic quantum numbers. Flavor nonsinglet hybrids 
with exotic quantum numbers are especially interesting because they cannot mix either with 
ordinary mesons or with glueballs. They can, however, mix with four quark (qqqq) states. 
In the quenched approximation mixing of hybrids and four quark states through sea quark 
loops (Fig. |I]a) is not present. However, the hybrid interpolating operators may still couple 
to four quark states through "hairpin" diagrams (Fig. |l|b). This coupling can be investigated 
by using four quark source operators (Fig. |^), and we find that it can even be a useful tool 
in computing the mass of the exotic mesons. 
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Figure 1: Quark line diagrams showing mixing with sea quarks (top), "hairpin" diagrams 
mixing hybrid and qqqq states in the quenched approximation (center), and an off diagonal 
propagator with a four-quark source and a hybrid sink (bottom). The vertical line indicates 
two meson operators at the same Euclidean time but different spatial coordinates. 
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Because of the partial suppression of mixing with four quark states and because we have 
a supply of quenched lattices at several lattice spacings, we have calculated exotic meson 
propagators in the quenched approximation at Q/g^ = 5.85 and 6.15. We have done this 
at a set of quark masses greater than or about equal to the strange quark mass, and at 
approximately the charmed quark mass. 



2 Hybrid operators and propagator computations 

To make an operator which creates a hybrid meson, we combine a quark, an antiquark and the 
color electric or magnetic field to form a color singlet with the desired spin, parity and charge 
conjugation. We construct these operators by combining representations of the continuum 
rotation group. An alternative approach using the symmetry group of the hypercubic lattice 
was presented by Mandulap, and is also developed in Ref. [^. 

Our hybrid operators have the generic structure tp'^Tip^F"'^ , where a and h are triplet 
color indices, F is some combination of Dirac matrices and derivatives, and F is the color 
electric or magnetic field, a color octet. Because we do not include "quark-line disconnected" 
diagrams in our propagator, all our meson propagators are flavor nonsinglets. 

The color electric and magnetic fields have J^^ = 1 and 1^^ respectively. The spin, 
parity and charge conjugation from the quark and antiquark are those of the available quark 
bilinears, listed here along with the corresponding mesons. The O"*"" bilinear, ip'^o'^j does not 
correspond to a qq state. Instead, it is the charge corresponding to a conserved current, the 
baryon number. Therefore we expect / d^x tpjoip \0) = 0. We may, however, calculate the 
propagator for this exotic operator, or use this bilinear as part of our toolkit for constructing 
hybrid operators. 
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We can also give the quark and antiquark a relative orbital angular momentum. This may 
be useful because in the nonrelativistic quark model the ai (1"^^), and hence the 0"*"" and 

hybrids constructed below, is a P wave state. The operator di= d i — d i inserted in the 
quark bilinear brings in quantum numbers 1 , where the negative charge conjugation comes 
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because C interchanges the quark and antiquark. Thus, a P-wave operator with ai quantum 

numbers, 1+"*", is eijkip'jj dk 'ip- This operator may be advantageous because it couples the 
"large" components of the quark spinor to the large components of the antiquark spinor. 



Figure 2: "pointlike" construction of F^jy. Each open loop represents the product of the 
links, minus the adjoint of the product. Each of these links may actually be a "smeared" 
link, as illustrated on the right side of the figure. 

For F^^ we use a "pointlike" construction, illustrated in Fig. |^. Each open loop represents 
the product of the links, minus the adjoint of the product. To improve the overlap of the 
operator with the hybrid meson, we can replace each link by a "smeared" link, as illustrated 
on the right side of the figure. The smeared link is the sum of the single link plus the 
three link paths displaced in the spatial directions, so there are four such staples for a link 
in a spatial direction and six for a link in the time direction. We have experimented with 
including staples displaced in the time direction in the smearing, and found that this distorts 
the propagators at short distances, and we will be including short distances in our fits. In 
our earlier work on two-flavor lattices we found that iterating the smearing twice gave a 
slightly better signal than a single smearing, but that four iterations of the smearing was 
worse. Therefore we have smeared twice in this work. We do expect that as the physical 
lattice spacing is decreased more iterations of smearing would be advantageous. 

Table |I| shows the various source and sink operators we have used. To construct a meson 
propagator, we first fixed the gauge to the lattice Coulomb gauge. We then constructed a 
wall source on one time slice, with a 1 at each lattice point for one color and spin compo- 
nent. A quark propagator was constructed by computing times this source. Then the 
wall source was multiplied by the source operator, which involved multiplication by Dirac 
matrices and components of the field strength. The result of this was used as a source for 
an inversion to compute the antiquark propagator. (This involves an extra 75 at each end 
of the propagator, from the standard identity = 'y^M*'-/^.) Finally, at each lattice point 
the antiquark propagator was multiplied by the desired sink operator, and dotted with the 
quark propagator. The result was summed over each time slice to get the zero momentum 
mesons. 



The first three operators in Table |1] are standard operators for the ^, 1 and l^"*" 
qq mesons. We will call these the "vr", "p" and "ai" respectively. However, our valence 
quarks are really much heavier than the physical u and d quarks, so they might be better 
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Table 1: Source and sink operators used for our propagators. The first column lists the 
name used for the operator in the text, and the second column lists the angular momen- 
tum, parity and charge conjugation. The third column is a shorthand indicating how the 
operator is constructed. In particular, the hybrid operators are constructed from one of the 
quark bilinear operators from the top block, combined with either the color electric or color 
magnetic field. The last column lists the actual operator. In this column, a and b are color 
indices, i, j and k Cartesian indices, and a, (3 and A are flavor indices included to indicate 
how the propagators are connected. 
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thought of as ss or, in one case, cc mesons. The fourth operator is a P-wave operator for 
the Oi- In the nonrelativistic quark model the ai is a P-wave state, so the operator with the 
spatial derivative will connect large components of the Dirac spinor to large components, 
and may give a better signal. This argument also applies to the hybrid operators that we 
will construct using this quark bilinear as a building block. 

The next set of operators are hybrid operators with the same quantum numbers as qq 
operators. In our previous work we verified, by computing propagators with a hybrid operator 
at one end and a qq operator at the other, that these hybrids mix with the corresponding 
qq operators (and don't mix with other quantum numbers) Q. The "mnemonic" column of 
the table indicates how the operator is constructed. For example, the first hybrid operator, 
which has pion quantum numbers, can be considered as a quark and antiquark in a 1 , or 
p state (but a color octet) combined with a color magnetic field, which has J^'" = V , to 
make a J = color singlet object. 

The remaining sections of Table || contain the operators with exotic quantum numbers. 
We have experimented with three 0^ operators. The first two are formed from the Oi quark 
bilinear and the color magnetic field, while the third is the Jb bilinear. Of these, the Op~ 
source and sink gave the best signal. (The Op~ operator in the form needed for the lattice 
computation is shown in the appendix.) 

There are two operators, using the "pointlike" and the "P-wave" ai bilinear respec- 
tively, and the color electric field. Again, we found that the P-wave operator gave a better 
signal. 



Finally, there are three 1 operators. The first is a quark and antiquark in a "p" state 



with the color magnetic field. As is well known [|10|, a 1 ^ hybrid can be constructed with 
the quark and antiquark in a relative S-wave state, and this is one of several arguments 
leading us to expect that it will be the lightest exotic hybrid. The second 1~+ operator is 
the charge bilinear ip^Qip combined with the color electric field. Of these two operators, we 
find that the first gives a better signal. The last operator is a four-quark operator. To use 
this operator as a source, we begin with the usual wall source, which is used as the source 
for a quark propagator computation. We then multiply the wall source by the pion operator 
(multiply by 75), and use this as the source for a computation of an antiquark propagator. 
The resulting propagator on the source time slice is then multiplied by the ai operator, 
(multiply by 7574), and the result is used as the source for another antiquark propagator. 
Obviously this is expensive, since it involves an extra propagator computation. We have 
not used this operator as a "sink" operator, since to do this we would need a separate extra 
inversion at each time slice, instead of just at the source time slice. However, the "crossed" 
propagator with the four quark operator as a source and the hybrid l""*" operator as a sink 
turns out to be useful. 
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3 Simulation results 



We used quenched lattices with the standard plaquette gauge action on 20^ X 48 lattices 
at Q/g"^ = 5.85 and on 32^ x 64 lattices at Q/g"^ = 6.15. These lattices were generated for 



quenched spectrum studies [pH 



At 6/g^ = 6.15 we evaluated propagators at five values of the Wilson hopping parameter, 
with the largest one chosen at approximately the charm quark mass. We used vr, p, oi, and 
all of the exotic operators in Table as source operators. For each source, we used all of 
the sink operators with the same quantum numbers, except for the operator. For each 
lattice, we used four source time slices. For the lightest quark mass exotic propagators, we 
used 30 lattices, with fewer lattices for the smaller k values. Because we do not need as 
many lattices to get good values for the vr and p masses, and because we did not implement 
the source until the project was already started, we do not have all the propagators on 
all of the lattices. 



At 6/(?^ = 5.85 we evaluated propagators with vr, p, ai(P), 1~+ and Q"^ sources. Since 
our previous work on G/g"^ = 5.6 two fiavor lattices and our concurrent work at Q/g^ = 6.15 
had found the 1~"*" to be the lightest of the exotics and the one for which we had the best 
signals, we did not do the 0"*"" and propagators on the 6/g^ = 5.85 lattices. We used 23 
lattices, with four source time slices on each lattice. 

Since propagators with different values of k were computed on the same quenched lattices, 
the mass estimates are strongly correlated. Also, there is the possibility of correlations 
among the propagators with different source time slices on the same lattice, which we have 
ignored in computing the covariance matrix. On the other hand, the different lattices are 
uncorrelated, unlike the situation in most full QCD calculations. To account for these 
correlations, especially when extrapolating or interpolating masses to different k, values, a 
jackknife analysis is useful. At Q/g'^ = 5.85 all of the propagators were run on the same set 
of lattices. There we used a single elimination jackknife, in which mass fits to the selected 
distance range and the extrapolation to Hc were repeated, each time omitting one lattice. 
For the mass estimates at Q/g"^ = 5.85 in Table § the first parenthesized error estimate 
is from the covariance matrix, and the second is from the jackknife analysis. These two 
estimates agree well. At Q/g'^ = 6.15 we do not have all the propagators or all the kappa 
values on every lattice. Therefore we divided the lattices up into five jackknife blocks, with 
each block containing about the same number of each kind of propagator, and ran the fits 
for the selected distance ranges, each time omitting one fifth of the lattices. Again, when 
there is a second parenthesized error estimate in Table ^ it is from the jackknife analysis. 

We first need an estimate of the lattice spacing. This can be done by extrapolating the 
p mass to the physical quark mass, or essentially to k^, or from the IS-IP mass splitting 
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for heavy quarks. Table ^ contains mass estimates for the pseudoscalar, vector and 1 + + 
qq ("vr", "p" and "cti") states. These mass estimates are crude by today's standards, but 
we need them only for approximately determining the lattice spacing. For = 5.85 we 
include vr and p mass estimates at additional k values, done on a 12^ spatial lattice, coming 
from our fs calculations |T^. 

At Q/g"^ = 5.85 we estimate I/kc = 6.205(3) from extrapolating the squared pion masses 
from the largest four k values. If we estimate the lattice spacing by extrapolating the p 
mass to Kc, we find = 1.80(6) GeV or 2.00(11) GeV, depending on whether nip or is 
extrapolated linearly in These extrapolations are plotted in Fig. ^ 

At 6/(7^ = 6.15 we do not have the luxury of extra pion masses near Kc- A linear fit of 
in 1/k does not work well, so we fit the pion mass at the four largest k's to a quadratic in 
1/k to estimate I/kc = 6.4895(10). This fit had a of 0.5 with one degree of freedom. An 
extrapolation of nip to this Kc gives amp{Kc) = 0.274(4), or = 2.81(4) GeV, where the 
error is a jackknife estimate. These extrapolations are plotted in Fig. ^. Another possibility, 
following the Fermilab group, is to use the splitting between the average S wave charmonium 
masses and the P wave masses to estimate the lattice spacing. At k, = 0.1350, using the l^'^ 
meson as the P wave mass, and |mi — + jniQ-+ as the S wave mass, with 457 MeV as the 
experimental value for charmonium, we get = 2.85(1) GeV. 

Fig. ^ shows exotic propagators at Q/g'^ = 6.15 and k, = 0.1350. For the 1 exotic we 
show three propagators. One has the 1 operator as its source and sink, the second the l""*" 
operator as the source and the 1^"*" operator as the sink, and the third has the operator 
as its source and the 1 operator as its sink. Fig. ^ is a similar figure for our smallest quark 
mass, K = 0.1520, except that instead of the 1^+ I2 ^ propagator we have a Iq^^ I2 ^ 
propagator. Compared to the propagators for conventional mesons, these exotic propagators 
are quite noisy, and we must use fairly short distances for the mass fits. As seen in Figs. ^ 
and ^, this problem becomes worse as the quark mass is made lighter. In fact, these exotic 
propagators fall below their statistical error at a distance smaller than the minimum distance 
we use for a single mass fit for the tt or p. This means that contamination by excited states 
might be a serious problem. We have therefore done fits to two exponentials for the exotic 
states. We can do this only over a small range of minimum distance. In particular, if we take 
the minimum distance too large we get fits with a very large excited state mass, which is 
essentially just a delta function removing the shortest distance point from the fit, and giving 
the same result as a one mass fit with a minimum distance one unit larger. Of course, we 
reject such fits since they are basically one mass fits. 

In these fits, we would like to use information from the different source and sink operators, 
by simultaneously fitting to two or more combinations of operators with different amplitudes 
for each source and sink but the same masses for all of them. 

{Oi{0)Oj{t)) = A^A^e-"^^' + A\A]e-'^'' . . . , (1) 
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Figure 3: Particle masses at 6/(7^ = 5.85, and extrapolations to The octagons are 

pion masses, and the line is a fit to linear in l/n. Squares and crosses are p and ai 
masses respectively. For the p and ai we show two extrapolations of the mass to Kc, one 
with mp^ai linear in I/k and the other with rrij^^^ linear in Finally, the bursts are the 
1 ^ exotic meson. We show two fits of the l""*" meson for each k. The lower one uses the two- 
source, two- mass fits as illustrated in Fig. ^ while the upper one uses only the l""*" 
propagator. 
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Figure 4: Particle masses at 6/(7^ = 6.15, and extrapolations to k^^. The symbols are the 
same as in Fig. ^, and the fancy plusses are a 0"*"" exotic. Again we show two fits of the 1""*" 
meson for each k, the lower one fitting both source operators simultaneously and the upper 
one fitting only the 1^^ 1 ^ propagator. The k = 0.1350, or I/k = 7.407, charmonium 
point is not in the range shown here. 
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Figure 5: Propagators for the l""*" exotic meson at Q/g"^ — 6.15 and k — 0.1350. The 
octagons are for a 1"+ source and sink, the diamonds for a l""*" source with a V^'^ sink, and 
the squares for a source with a 1~+ sink. 
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Figure 6: Propagators for the 1""*" exotic meson at &/g'^ = 6.15 and k = 0.1520. The 
octagons are for a l""*" source and sink, the diamonds for a I2''' source and sink, and the 
squares for a Q'^ source with a sink. Plusses indicate the absolute value of a propagator 
that changes sign with increasing distance. 
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Figure 7: Effective masses for the p and 1~+ exotic at &/g^ — 6.15 and k — 0.1350, 
approximately the charm quark mass. The octagons are the effective mass from propagators 
with the 1^+ operator as both source and sink, and the squares are the effective mass with 
the Q'^ source and l""*" sink. These propagators were fit with two source amplitudes and two 
masses as described in the text, over the distance range 2 to 11. The sohd hues near the plot 
symbols are the effective masses reconstructed from the fit. The upper pair of horizontal 
lines indicates the ±la range for the ground state mass in this fit. The diamonds are the p 
(more accurately, the ip) effective mass. The horizontal bars near the diamonds are the ilo" 
limits on the ground state mass from fits to the p propagator. The lines running from d = 2 
to 6 are a two mass fit to the propagator over this range, while the closely spaced lines for 
d > 6 are from a single exponential fit over distance range 21 to 25. 
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Figure 8: Effective masses for the 1 meson a.t = 6.15 and k = 0.1520, approximately 
the strange quark mass. The symbols and lines have the same meaning as in Fig. 0. These 
effective masses correspond to the propagators in Fig. |[ 
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Figure 9: Effective masses for the 1 meson and the p at G/g"^ 
The symbols and hnes have the same meaning as in Fig. 0. 



5.85 and k = 0.1450. 
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Figure 10: Effective masses for the 1 meson and the p at Q/g"^ 
The symbols and hnes have the same meaning as in Fig. 0. 
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Table 2: Mass estimates for ordinary qq mesons. The Oi mass estimates used the ai(P) 
source and sink. An asterisk indicates a point from a 12^ lattice. Where present, a second 
parenthesized error is a jackknife estimate. Fits at «; = 0.154 for Q/g^ = 6.15 are jackknife 
extrapolations to Kc- 
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1 


3-9 


3.9/5 


1.81(3) 






1-+ ^ 


1-+ 


1 


4-10 


3.5/5 


1.88(8) 






Q'- 


1-+ 


1 


3-7 


0.7/3 


1.65(5) 








1-+ 


2 


1-9 


7.3/5 


1.72(9) 


2.45(19) 




Q^l-+^ 


1-+ 


2 


1-8 


9.1/10 


1.71(3) 


2.42(6) 






1-+ 


2 


2 8 


6.6/8 


1.75(1)(5) 


2.86(44) 


0.1500 


1-+ - 


1-+ 


2 


1-9 


5.8/5 


1.54(8) 


2.33(19) 






1-+ 


2 


1-8 


7.4/10 


1.47(3) 


2.20(4) 






1-+ 


2 


2-8 


6.2/8 


1.51(5)(6) 


2.31(21) 


0.1525 


1-+ - 


1-+ 


2 


1-9 


5.2/5 


1.45(8) 


2.28(18) 




Q^ 1-+ ^ 


1-+ 


2 


1-8 


7.4/10 


1.33(3) 


2.10(4) 




Q^ 1-+ 


1-+ 


2 


2-8 


6.0/8 


1.38(6)(7) 


2.08(15) 


0.1550 


1-+ - 


1-+ 


2 


1-9 


5.0/5 


1.36(9) 


2.22(17) 




gM-+- 


1-+ 


2 


1-8 


9.3/10 


1.19(4) 


2.00(3) 






1-+ 


2 


2-8 


6.0/8 


1.20(7)(8) 


1.88(10) 



Table 3: Mass estimates for the exotic 1""*" meson for = 5.85. Where two source 
operators are hsted, a simultaneous fit was done to propagators from both sources, with 
the masses forced to be the same for each source. Where two masses were used in the fit, 
the last column shows the excited state mass produced by the fit. Where present, a second 
parenthesized error is a jackknife estimate. 
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K 


Source(s)- 


-i>Sink 


masses 


Fit Range 


XVDOF 


Mass 


M* 


0.1350 


1"+ - 


1"+ 


2 


1-11 


2.3/7 


1.61(2) 


2.52(4) 




1"+ - 


1"+ 


2 


2-11 


1.3/6 


1.58(3) 


2.42(11) 




1"+ - 


1"+ 


2 


3-11 


1.2/5 


1.59(4) 


2.54(42) 






1 + 


2 


1-11 


18.1/16 


1.63(1) 


2.57(3) 




Q\l-+ - 


1"+ 


2 


2-11 


6.8/14 


1.58(3) 


2.38(9) 




Q\l-+ - 


1"+ 


2 


3-11 


6.5/12 


1.57(4) 


2.28(23) 


0.1450 


-I _L 

1"+ - 


1"+ 


2 


1-11 


11.5/7 


1.189(16) 


2.19(4) 




1"+ - 


1"+ 


2 


2-11 


4.5/6 


1.143(24) 


2.24(7) 






1 + 


2 


2-11 


15.2/14 


1.109(16) 


1.90(4) 




Q^l +- 


1 + 


2 


3-11 


8.8/12 


1.134(19)(34) 


2.04(18) 






1 + 


2 


4-11 


4.2/10 


1.161(39) 


1.62(47) 


0.1480 


1"+ - 


1"+ 


2 


2-11 


4.8/6 


1.013(22) 


2.16(6) 




1"+ - 


1"+ 


2 


3-11 


3.0/5 


1.039(25) 


2.30(43) 




^^A -1 _l_ 


1 + 


2 


2-11 


30.4/14 


0.955(15) 


1.76(4) 




/04 1 h 


1 + 


2 


3-11 


"1 r' 1 /to 

15.1/12 


0.980(21)(36) 


1.74(12) 




QM-+- 


1-+ 


2 


4-11 


6.3/10 


1.019(29) 


1.67(33) 


0.1500 


1-+ - 


1-+ 


2 


2-11 


8.8/6 


0.937(24) 


2.16(6) 




1-+ - 


1-+ 


2 


3-11 


8.7/5 


0.943(29) 


1.91(24) 






1-+ 


2 


2-11 


48.9/14 


0.868(17) 


1.70(03) 




QM-+- 


1-+ 


2 


3-11 


22.7/12 


0.897(22)(19) 


1.63(10) 




QM-+- 


1-+ 


2 


4-11 


15.5/10 


0.943(26) 


1.88(42) 


0.1520 


1-+ - 


1-+ 


2 


2-11 


6.9/6 


0.839(26) 


2.13(6) 




1-+ - 


1-+ 


2 


3-11 


6.4/5 


0.820(37) 


1.63(16) 






1-+ 


2 


2-11 


38.7/14 


0.789(23) 


1.66(4) 






1-+ 


2 


3-11 


16.0/12 


0.796(33) (47) 


1.54(11) 




QM-+- 


1-+ 


2 


4-11 


14.3/10 


0.827(40) 


1.76(50) 


0.154 


QM-+- 


1-+ 


2 


extrap. 




0.705(na)(32) 





Table 4: Mass estimates for the exotic 1 meson for 6/(?^ = 6.15. The format is the same 
as Table ^. The final line is an extrapolation to using the distance range 3-11 for the 
four largest k values. 
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K 



Source(s)^Sink masses Fit Range xV^OF Mass 



0.1350 



0.1450 



0.1480 



0.1500 



0.1520 



0+- 
0+- 
0+- 

0+- 
0+- 

op 
0+- 

ot 
0+- 
0+- 

w 



0+- 
0+- 
0+- 
0+- 



oy- 

0+- 
0+- 
0+- 

"or 

oj- 

0+- 
_oJi 
0+- 

or 

0+- 
0+- 

0+- 



0+- 
0+- 
0+- 
0+- 



1 
1 
1 

_2_ 

T 
1 
1 

_2_ 

T 
1 
1 

2 

T 
1 
1 

_2_ 

T 
1 
1 

2 



5- 9 

6- 10 

7- 11 
2-11 

"5^ 

6- 10 

7- 11 
2-11 

5- 9 

6- 10 

7- 11 
2-11 

5- 9 

6- 10 

7- 11 
2-11 

5- 9 

6- 10 

7- 11 
2-11 



3.6/3 
2.7/3 
1.2/3 
4.5/6 
2.9/3 
3.6/3 
2.7/3 
3.9/6 
1.9/3 
2.2/3 
1.9/3 
3.5/6 
0.7/3 
5.4/3 
4.9/3 
7.5/6 
2.7/3 
1.9/3 
1.7/3 
4.2/6 



1.72( 

1.69( 

1.63( 

1.66( 

1.32( 

1.30( 

1.23( 

1.27( 

1.20( 

1.18( 

l.lll 

1.16( 

1.131 

l.OSi 

1.01 

l.lOl 

1.08( 

0.98( 

0.94( 

1.04( 



M* 



(5) 2.52(6) 



(5) 2.24(7) 



)(3) 2.16(6) 



2] 
4) 

9' 
3) 

3) 
6) 
13) 
3)(3) 



(2) 2.16(6) 



2.13(7) 



Table 5: Mass estimates for the exotic O"*" meson for G/g"^ = 6.15. The format is the same 
as Table 
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where i and j label the source and sink operators respectively, and mo and mi are the ground 
state and excited state masses. For this to be useful, the relative overlaps of the different 
operators with the ground state and excited state should be as different as possible. For 
the 1""*" propagators it turns out that the 1 and l^"*" operators have essentially the same 
effective masses. Since the propagators with l^"*" operators are noisier than those with l""*" 
source and sink, including these propagators in the fitting did not help (any gain in statistics 
was not worth the extra degrees of freedom in the fitting.) However, the correlator generated 
from the Q'^ source, which we introduced to investigate coupling to four quark states, does 
have an effective mass at short distance significantly different from that generated from the 
hybrid 1 source. This can be seen in Fig. 0, which shows the effective masses corresponding 
to the 1^+ 1^+ and Q'^ — > 1^+ propagators at our heaviest quark mass, k = 0.1350, and in 
Fig. ^ which shows the effective masses at our lightest quark mass, k, = 0.1520, corresponding 
to the propagators in Fig. |^. As might be expected from this, simultaneously fitting the l""*" 
and source propagators, each with the 1""*" sink operator, to two masses gave the best 
mass estimates. The ground state mass from this fit, and the effective masses from the fit, 
are also shown in Figs. |^ and ^. Of course, this two source fit makes the assumption that 
the excited state (or combined effect of many excited states treated as a single state in the 
fitting program) is the same in both propagators. Therefore we tabulate results both from the 
simultaneous fits to two source operators and from fits using only the 1 source operator. 
In the cases of the O"*"" and propagators, we have not investigated four quark source 
operators. For the 0"*"~ and hybrid operators, the P wave source and sink operators 
generally gave the best statistical errors, and so were the only ones we fit. 

One might still worry that we are not extracting the correct ground state masses from such 
short distances. As a partial check, we take the p propagator, for which a quite convincing 
plateau in the effective mass is seen at larger distances, and make a two mass fit to this 
propagator at the same distances we use for the hybrid fits. (Since the hybrid fits are 
dominated by the more accurate points at the smaller distances in the fitting range, we 
fit the rho mass to the smaller distances in the range.) In this case, at Q/g"^ = 6.15 and 
K = 0.1350, which is the charm quark mass, we find that the p (more accurately, the J/ip) 
mass from a single exponential fit to distance range 21 to 25 is 1.1067(5) {x^/d = 2.4/3), 
while a two mass fit over distance range 2 to 6 gives a ground state mass 1.109(8) with 
an excited state mass of 1.264(21) {x^/d = 1.5/1), in excellent agreement with the single 
exponential fit from long distances. The p (J/V') effective mass and these fits are also shown 
in Fig. ^ While this result is encouraging, we should caution the reader that the number 
of excited states and the mass gap between the ground and excited states might be very 
different for the p meson and the exotic mesons. 

Another important test is to verify that the mass estimates are independent of the fitting 
range used. Here we are not in as good a position, since we generally have only two or three 
minimum distances where we can get a two mass fit with reasonable x^- However, within 
the fairly poor statistical errors, the exotic mass estimates are generally consistent among 
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these fits. 



The 6/(7^ = 5.85 fits to the 1 are done in similar fashion. In Figs. ^ and [1^ we show 
the effective masses and the fits to the propagators for our largest and smallest k values, 
n = 0.1450 and 0.1550, at = 5.85. 

Tables |^ and § contain selected mass fits for the 1 mesons at Q/g'^ = 5.85 and 6.15 
respectively. When two masses were used in the fit, both the ground state and the excited 
state mass are tabulated. However, this excited state mass is almost certainly some sort of 
weighted average of many states, and should not be taken seriously as a mass estimate. 

Finally, Table ^ contains selected fits for the 0^ exotic at Q/g"^ = 6.15 This particle is 
clearly heavier than the l""*", and our estimates for its mass are worse. This is partly because 
we have only the one source operator for this meson, and partly because (nonrelativistically) 
the quark and antiquark are in a relative P wave state, and we were led to use a more 
complicated source operator. 

We were unable to get credible mass estimates from the propagators, suggesting that 
this state, if it exists at all, is even heavier than the 1 ^ and 0^ . 



4 Discussion and Conclusions 



In this concluding section we discuss the conversion of our mass estimates from lattice units 
to physical mass units, estimate systematic errors, and describe briefly the observational 
status of 1~+ exotic hadrons. 

Figures |^ and ^ collect and display our results for nonexotic mesons and for the l""*" 
and the O"*"" exotic mesons. We also show extrapolations to Kc, where the error on the 
extrapolations comes from the jackknife analysis. For the 1 hybrid we plot fits to both the 
one source operator and two source operator mass estimates. The difference between these 
fits is an indication of the possible systematic error from excited states in the propagators. 
While this difference is small for the charmonium point, at the strange quark mass for 
Q/g'^ = 6.15 it amounts to 117 MeV, and when extrapolated to Kc becomes 165 MeV. 

One of the known problems with the Wilson quark action is that it consistently under- 
estimates spin splittings of hadrons[|13|. This suggests that the average mass of the S wave 
mesons might be a better mass standard than the p mass alone. (We have already used this 
logic in section 3 when we used the splitting between the P-wave and the average S-wave 
charmonium mass as a length scale.) In Fig. we plot the ratio of the 1"'"+ P-wave meson 
(oi) to the average S-wave meson mass, |mi — -|- |mo-+. For the horizontal scale we use 
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2.5 




1 2 3 4 5 

ms/(0.75*m^+0.25*m„) 



Figure 11: Ratios of the 1"'""'" (P wave) meson mass and the 1""*" exotic mass to the average 
S wave meson mass. Diamonds and bursts are the 1^^ at 6/(7^ = 5.85 and 6.15 respectively, 
and squares and octagons are the 1^+ at 5.85 and 6.15 respectively. The horizontal and 
vertical scales and the vertical and horizontal lines are described in the text. 
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the average S-wave mass divided by the average S-wave mass at the physical point, where 
mT^/nip = 0.18. In other words, the units of the horizontal axis are |mp + ^m^ = 610 MeV. 
The vertical scale is the ratio of the P-wave l"*""*" meson mass or the l""*" exotic meson mass 
to the S-wave mass. The vertical lines indicate the strange quark and charmonium points. 
In locating the strange quark line, we have followed the UKQCD procedure of using an "un- 
mixed" rjstrange uiass of 680 MeV0. By definition, the light quark point is the left hand side 
of the graph. Where they intersect the left side of the graph and the two vertical lines, the 
three bold horizontal lines indicate the experimental values of the l^"*" meson mass divided 
by the S-wave mass for light quarks, strange quarks and charm quarks respectively. 

A nice feature of this graph is the good agreement of the P-wave masses between the 
G/g"^ = 5.85 and 6.15 lattices. For heavy quarks the agreement of the l""*" exotic is equally 
good, but there is some difference for the light quarks. We suspect that the Q/g"^ = 5.85 
points are incorrect here. In Fig. we see that the l""*" fit at Q/g^ = 5.85 and k = 0.1550 
is questionable at best, and Fig. ^ or Table |^ show that had we used the one source fits for 
this point we would have obtained a 13% larger mass. The exact agreement between our 
predicted and the observed charmonium mass is undoubtedly fortuitous, but this agreement 
and the trend toward agreement between the 5.85 and 6.15 estimates at larger quark mass 
encourages us to quote a charmed hybrid meson mass with suitable caveats. 



Using this average S-wave mass as the length scale for the charmonium exotics gives 
an uncorrected mass of 4390(80) MeV for the 1"+ and 4610(110) MeV for the 0+". We 
emphasize that these quoted errors are statistical only, and do not take into account con- 
tamination from excited states in the chosen fitting ranges, or discretization errors in the 
gauge and quark action, or effects of the quenched approximation. As mentioned above, 
the errors from excited state contamination, which can be crudely estimated by looking at 
how the estimated mass varies as a function of minimum distance in the fit or whether the 
second source operator is included, are probably about the same size as the statistical errors. 



The close agreement of the l^"*" masses and the 1 masses for heavier quarks in Figure 11 



suggests that errors from the nonzero lattice spacing might be small. However, we also note 
that in our studies of charmed pseudoscalar meson decay constants we find a discretization 
error in the decay constants of 10% to 15% at G/g"^ = 6.15|[T^. It might be that decay 
constants are more sensitive to lattice spacing errors than mass, since the decay constants 
are basically wave functions at the origin, and are strongly affected by the coarseness of the 
lattice at short distance. Although it is little better than a quess, we propose using 15% 
of the splitting between the hybrid state and the cc states, or 200 MeV, as an estimate of 
the systematic error from both excited state contamination and lattice artifacts. The largest 
systematic error probably comes from using the quenched approximation, and this is the 
hardest error to estimate. It will surely be large, since, as discussed in the introduction, 
these hybrid states can mix with four quark states. In this work we looked at mixing with 
states containing four heavy quarks, but in the real world the important four quark states 
would contain the charmed quark and antiquark and a light sea quark and antiquark. This 



25 



charmonium 1 is quite far above the DD threshold. Because of the remaining systematic 
uncertainties, it is not clear whether it is above the S wave + P wave DD threshold, which 



in many model calculations determines whether its decay width is large] 15 



At the strange quark mass (our largest k at Q/g"^ = 6.15), we estimate the mass of the 
1^^ hybrid to be 2170 ± 80 ± systematic MeV. As noted at the beginning of this section, the 
systematic error from excited state contamination is on the order of 100 MeV, and the error 
from nonzero lattice spacing is probably as large or larger. Considering the large errors, 
this estimate is consistent with the mass quoted by the UKQCD collaboration, 2000(200) 
MeV0. 

If we take seriously the extrapolation of the l""*" mass to light quarks in Fig. ^ we get a 
mass of 1970(90) MeV for the light quark exotic hybrid, again with large systematic errors. 
Considering the lack of agreement in the l""*" mass estimates at Q/g"^ = 5.85 and 6.15 at light 
quark mass and the effects of extrapolating in kappa, we use 300 MeV as an estimate of the 
error on this number from lattice artifacts. 



We briefly review the observational evidence relevant to our results. The particle data 
table does not list any confirmed 1 hybrid meson states. However a number of 
potential candidates are mentioned. There is some evidence for a 1~^ hybrid state with a 
mass of around 1.4 GeV, and another with a mass around 1.9 GeV. The original evidence 
for a 1^^ state at 1.4 GeV found by the GAMES collaboration [|l^ was criticized in [|18|. 
However recent work by the E852 collaboration reports evidence for a l""*" hybrid state 
with a mass of 1370 ± IGtll MeV. This paper Ijl^ also lists other experiments that have 
reported a low mass for the 1 state. Our result favors a hybrid assignment for states 
around 1.9 GeV 
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However we stress that more simulations are required to quantify 
and reduce the systematic errors in our results, before definitive results for the mass of the 
1 state can be obtained from quenched lattice QCD. 



In Table ^ (obtained from ) we collect results for the mass of the light 1^^ hybrid 
obtained from a variety of models. We note that our prediction for the mass of the l""*" state 
is consistent with the flux tube estimate. Only the bag model calculation obtains a mass 
close to the E852 experimental result. 

Unfortunately, no good candidates for the theoretically expected 1 ^ hybrid charmonium 
state have been observed [1^. Such a state is also expected in flux tube models, which are 
most plausible for heavy quark systems and furthermore provide information about branching 
ratios. Table ^also contains some predictions (taken from 0) for the mass of the 1 state 
in the charmonium system. 

Although we have concentrated on exotic hybrids in this work, lattice methods, partic- 
ularly with dynamical quark loops included, could also address interesting issues regarding 
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Mass (GeV) 



Method 



Light quark 1 + mass 
1.3 ^ 1.8 
1.8 ^ 2.0 
1.8 ^ 1.9 
2.1 2.5 

quench ) 



Bag model 

Flux tube model 

Flux tube model of Close et 

QCD sum rules (mostly after 1984) 

This work 



Charmonium 1 mass 

^ 3.9 

4.2 ^ 4.5 

4.1 ^ 4.2 

4.19 ± syst. 

4.1 ^ 5.3 

4:.39{8stat.) {'iOlattice) { • • quench ) 



Adiabatic bag model 

Flux tube model 

Flux tube model of Barnes et al. 

Heavy quark LGT§ 

QCD sum rules (mostly after 1984) 

This work 



Table 6: Predictions for the mass of the light and charmonium 1 + states from various 
approaches to QCD, obtained from |]^. 



nonexotic states. For example, Close and Page characterize the '?/'(4040) and ?/'(4160) states 
as mixtures of cc and cc5'[^. Close proposes that hybrid states could help explain the 
"anomalous" production of charmonium observed by the CDF group |2^. In both cases the 
occurrence of a nonexotic hybrid state in the range 4 to 4.3 GeV is essential. 



Future efforts in lattice calculations of exotic hybrids should be directed at reducing the 
variety of systematic errors, including finite lattice spacing, excited state contamination, and 
the effects of quenching. Particularly important is to understand the extent of mixing with 
four-quark states. 
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Appendix 

Here we write the 0"^~ P-wave operator, Op", in the form needed for the lattice computation. 

We are interested in zero spatial momentum, so wc arc summing over spatial coordinates, 
and can freely translate the summation variable. Wc wish to write the operator as iplx) 
times a sum of fields at x and neighboring points. In these expressions, it is understood that 
fields at neighboring points must be parallel transported to x, so that 'ip{x)ip{x + k) means 
'ip{x)Uk{x)ip{x + k). 



X jk 

= EE i>\x)^jr{s+k)Ff,{x) 

X jk 

- i,\x)^,r{x-k)Ff,{x) 

- i,\x + kYi,r{x)Ff,{x) 
+ i,\x-k)^,r{x)Ff^{x) 

= EE i^\x)ijr{x+k)F${x) 

X jk 



- ^\x)^,r{x-k)Ffi:{x) 

- ^\x)^^r{x-k)Ff^{x-k) 
+ i)\x)^jij\x + k)F^t{^ + k) 

E E ^'(^)7. {r{x + k) {f${x) + F${x + k))) 



X jk 



k) {F^^ix) + F^^ix - k))) (2) 



Similarly, for the P-wave source: 

= E E ^^Jk i^S i^k - ^k) rFst 

X ijk 

= EEe^.fc( i^\xhjrix+ k)F^Kx) 

X ijk 

- i,\x)^,rix-k)F,^^ix) 
28 



X ijk 



X ijk 



- i,\x)^,r{x-k)FSt{x-k) 
+ i)\x)-i,ij\x + k)FS^{x + k) 



- (^'^(x-^)(Fof(f)+Fof(f-fc) 



(3) 
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